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ABSTRACT 

Gamma-Ray Bursts (GRBs) are widely thought to originate from collimated jets of material moving 
at relativistic velocities. Emission from such a jet should be visible even when viewed from outside 
the angle of collimation. I summarize recent work on the special relativistic transformation of the 
burst quantities Aj so and A pea k as a function of viewing angle, where Aj so is the isotropic-equivalent 
energy of the burst and Apeak is the peak of the burst spectrum in the power vF^. The formulae 
resulting from this work serve as input for a Monte Carlo population synthesis method, with which 
I investigate the importance of off-jet relativistic kinematics as an explanation for a class of GRBs 
termed “X-ray Flashes” (XRFs). I do this in the context of several top-hat shaped variable opening- 
angle jet models. I find that such models predict a large population of off-jet bursts that are observable 
and that lie away from the A pea k oc A iso relation. The predicted burst populations are not seen in 
current datasets. I investigate the effect of the bulk 7 value upon the properties of this population of 
off-jet bursts, as well as the effect of including an n 0 -A 7 correlation to jointly fit the A pea k oc E-J o and 
A P eak oc A^ relations, where Sl 0 is the opening solid angle of the GRB jet. I find that the XRFs seen 
by HETE-2 and BeppoSAX cannot be easily explained as classical GRBs viewed off-jet. I also find 
that an inverse correlation between 7 and flo has the effect of greatly reducing the visibility of off-jet 
events. Therefore, unless 7 > 300 for all bursts or unless there is a strong inverse correlation between 
7 and H 0 , top-hat variable opening-angle jet models produce a significant population of bursts away 
from the A pea k oc Aj 1 / 2 and A pea k oc Al® relations, in contradiction of current observations. 

Subject headings: Gamma Rays: Bursts — ISM: Jets and Outflows — Shock Waves 


1. INTRODUCTION 

The importance of collimated jets in GRBs was high¬ 
lighted by the extremely large isotropic-equivalent en¬ 
ergies (Ai so ) of very luminous events like GRB 971214 
(Kulkarni et al. 1998) and GRB 990123 (Kulkarni et al. 
1999) and by the observation of breaks in afterglow light- 
curves (Rhoads 1997; Sari, Piran & Halpern 1999; Har¬ 
rison et al. 1999). Frail et al. (2001) and Bloom et al. 
(2003) corrected the isotropic-equivalent energies by the 
beaming fraction obtained from the jet break time in af¬ 
terglow light-curves and found that the values of the en¬ 
ergy released in 7-rays (A 7 ) were tightly clustered around 
10 51 ergs. Recently Ghirlanda et al. (2004) have shown 
that a tight correlation exists between A 7 and the peak 
of the vF v spectrum in the rest-frame, A pea k. Recent re¬ 
sults by Sakamoto et al. (2005a) obtained from HETE-2 
(Ricker et al. 2001) observations have shown that XRFs 
(Heise et al. 2001; Kippen et al. 2001), X-ray-rich GRBs 
and GRBs lie along a continuum of properties and that 

XRFs with known redshift extend the A pea k cx A iso rela¬ 
tion predicted by Lloyd-Ronning et al. (2000) and found 
by Amati et al. (2002) to over five orders of magnitude 
in A; so (Lamb et al. 2005). 

Relativistic kinematics implies that even a “top-hat”- 
shaped jet will be visible when viewed outside its angle 
of collimation; i.e., off-jet (Ioka & Nakamura 2001). Ya- 
rnazaki, Ioka & Nakamura (2002, 2003) used this fact to 
construct a model where XRFs are simply classical GRBs 
viewed at an angle 9 V > 9q, where 9q is the half-opening 
angle of the jet and 9 V is the angle between the axis of the 
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jet and the line-of-sight. The authors showed that such 
a model could reproduce many of the observed charac¬ 
teristics of XRFs. Yamazaki, Ioka & Nakamura (2004a) 
showed that in such a model, the distribution of both on- 
and off-jet observed bursts was roughly consistent with 
the A pea k cx E-J o relation. 

In this paper, I use the population synthesis method 
developed by Lamb, Donaghy, & Graziani (2005) and in¬ 
corporate the relativistic emission profiles calculated by 
Graziani et al. (2005), to predict the global properties of 
bursts localized by HETE-2 and BeppoSAX. I consider 
the possibility that the XRFs observed by HETE-2 and 
BeppoSAX are primarily regular GRBs observed off-jet 
(Yamazaki, Ioka & Nakamura 2004a) and show that it is 
difficult to account for the observed properties of XRFs 
in this model. However, since the effect of special rela¬ 
tivity on off-jet emission must exist, I seek to understand 
its relative importance in the context of current models 
of GRB jets. I revisit the top-hat variable opening-angle 
(THVOA) jet model put forward in Lamb, Donaghy, & 
Graziani (2005), now including the effects of relativis¬ 
tic kinematics on off-jet emission. I present results for 
several models which explore various regions of the pa¬ 
rameter space in 7, A 7 and 9$. 

For this paper, I only consider the effect of relativistic 
kinematics on off-jet emission 1 from uniform or “top-hat” 
jets; we will consider the effects on Fisher-shaped (Don¬ 
aghy, Lamb, & Graziani 2005a) and Gaussian-shaped 

1 In this paper I use the terms “off-jet emission” or “off-jet rela¬ 
tivistic kinematics”, rather than “off-axis beaming”, to emphasize 
that such emission is a direct consequence of special relativity and 
that it is primarily important beyond the edge of the jet. 
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(Zhang et al. 2004) jets in a future publication (Donaghy, 
Lamb, & Graziani 2005b). I describe my population syn¬ 
thesis method in §2 and present the results for various 
models in §3. I discuss the results in §4 and draw some 
conclusions in §5. Preliminary results were reported in 
Donaghy (2005a). 


2. METHOD 

2 .1. Off-Jet Relativistic Kinematic Formulae 

Relativistic kinematics causes frequencies in the rest 
frame of the jet to appear Doppler shifted by a factor 
S = 7(1 — /3 cos 9), where (3 is the bulk velocity of the jet 
and 8 is the angle between the direction of motion and 
the source frame observer. The simulations I describe 
chiefly deal with the kinematic transformation of two im¬ 
portant burst quantities, £) so and Ap ea k> as a function 
of viewing angle, 8 V . In the simplest “toy” model, these 
quantities transform as A pea k oc d _1 and A; so <x <5 , 

from which arises the relation A pea k oc A js ' o (Yamazaki, 
Ioka & Nakamura 2002). In the more complete model 
of Graziani et al. (2005), this relation is satisfied only in 
the limit 9 V 9q. 

The complete relativistic kinematic expressions involve 
convolution of the Doppler function and the intrinsic pro¬ 
file of the jet. For an arbitrary smooth profile, an efficient 
algorithm exists to calculate the profiles; for the case I am 
interested in, the uniform or “top-hat” profile, a closed 
analytic expression can be given. The formulae are de¬ 
rived in Graziani et al. (2005), and I summarize them be¬ 
low. The current model differs slightly from that used by 
Yamazaki, Ioka & Nakamura (2004a) in that we consider 
steady-state emission, rather than the evolution of burst 
properties due to time-of-flight effects. The model there¬ 
fore applies to burst-averaged data products like -E) so and 

A pea k- 

The observed isotropic-equivalent energy, A; so , of the 
jet as a function of 9 V is given by, 

jfj true 

Eiso = OR 4 TT -7TT [f(P - cos 9 V ) - f ({3 cos 9 0 - cos6» v )] 

zp 7 ^(l — cost/o) 

(1) 

where 

. _ 7 2 ( 2 y 2 - 1 )z 3 + ( 37 2 sin 2 9 v )z + 2 cosd v sin 2 9 V 
f{Z) = (z 2 + 7" 2 sin 2 9 V ) 3 / 2 ’ 

( 2 ) 

and A* rue is the total energy emitted by the jet in gamma 
rays and serves as the energy scale for the emission pro¬ 
file. 

The transformation of Speak as a function of 9 V is 
slightly more complicated. The detailed physics under¬ 
lying the prompt emission of gamma-ray bursts is not 
yet well understood. In particular, an explanation of 
the prompt emission spectrum (apparently universally 
parametrized by the Band function (Band et al. 1993)) is 
currently lacking. The observed spectrum (including the 
value of Speak) is almost certainly due to superpositions 
of emission from different regions on the jet, convolved 
with relativistic kinematic effects. A detailed explana¬ 
tion of how this forms a Band spectrum is beyond the 
scope of this paper. Instead we calculate the average 
Doppler shift across the jet as a proxy for A pea k. The 


average shift, (D), is given by, 

= -1 f{/3 - cos 6 > v ) - f{f3 cos 6< 0 - cos 9 V ) 
^ g{(3 — cosf? v ) — g(/3cos9o — cos9 v ) ’ 


where 


27 2 z + 2 cos 6 * v 
(z 2 + 7 -2 sin 2 9 V ) 1 / 2 


(4) 


Epeak is related to ( D ) via, 

E peak = (D)xE^\ (5) 

where E^^ is the (unknown) peak energy of the burst 
spectrum in the rest frame. We remove this unknown 
normalization by requiring all bursts to obey the Apeak oc 
1 /2 

EJ Q relation at the center of the jet, thereby fixing this 
normalization to that of Ei so . Thus, 


Apeak 


{DM -(on) 

(A)(0)) peak 


{DM 

{Dm 


(Ei° 0 n) /E A y\ 

( 6 ) 


where Ajg™ is simply Equation 1 evaluated at 9 V = 0. 
Figure 1 shows Aj so and (£>) plotted as functions of 9 V 
for various values of 9p and 7 , and Figure 2 shows the 
corresponding trajectories in the [Aj so ,(A)]-plane as 9 V 
increases away from the jet axis. 


2 . 2 . Monte Carlo Simulations 

The population synthesis Monte Carlo simulations de¬ 
scribed in this paper follow the method presented in 
Lamb, Donaghy, & Graziani (2005). Beginning in the 
rest-frame of each burst I specify 9q, A* rue and 7 by draw¬ 
ing from various input distributions (depending on the 
model used). A value for 9 V is drawn from the distribu¬ 
tion dN = sin(0 v ) d9 v . I then calculate Aj so and A pea k 
from Equations 1 and 6 . 

I also introduce two Gaussian smearing functions to 
add a stochastic element to the simulations. I draw val¬ 
ues for the coefficient C A in Equation 6 from a narrow 
lognormal distribution to model the observed width of 
^ 1 /2 

the Ap ea k oc A is ' 0 relation. Another narrow lognormal 
distribution is used to generate a timescale for each burst 
that converts fluences to peak fluxes. For both of these 
Gaussians I use the same parameters as in Lamb, Don¬ 
aghy, & Graziani (2005). 

I then draw redshift values from a model for the star- 
formation rate (Rowan-Robinson 2001), transform the 
burst quantities to the observer-frame and construct a 
Band spectrum (assuming a = —1 and f3 = —2.5). Using 
the observer-frame Band spectrum I can calculate photon 
and energy fluences and peak fluxes in any desired pass- 
band. By comparing with the peak photon flux thresh¬ 
olds as described by Band (2003), I determine if a burst 
would be detected by a given instrument. In this work, I 
primarily employ the detector thresholds from the WXM 
on HETE-2, scaled to include triggers on timescales up 
to 5 sec. 

One change in the method from Lamb, Donaghy, & 
Graziani (2005) is necessarily the treatment of off-jet 
events. In that paper, simulation of all off-jet events 
was bypassed by drawing from a power-law distribution 
in fi 0 with index (5 s ; m = <5t rU e — 1 ; as discussed in that 
work, this is equivalent to working in the 7 —» 00 limit. 
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Fig. 1.— Emission profiles (top row) and average Doppler shift (bottom row) as a function of viewing angle, for a range of 7 and Oo = 0.1 
rad (left) and 9g = 0.01 rad (right). Prom Graziani et al. (2005). 




Fig. 2. — Trajectories in the [i5; so ,(Z))]-plane (“Amati plot”), moving from right to left as 6 V increases, for do = 0.1 rad (left) and 
@0 = 0.01 rad (right). The axes are scaled so that a 45° line corresponds to a slope of 1/2. From Graziani et al. (2005). 


In that paper, we concluded that the data requested a 
model with approximately equal numbers of bursts per 
decade in all observed burst quantities, corresponding 
to (5 s i m = 1. Equivalently, in this paper I draw from a 
power-law distribution in fio with index Jtrue = 2, and 
simulate the distribution of viewing angles, 8 V , to deter¬ 
mine which off-jet events are detected. 

The vast majority of jets are observed from large view¬ 
ing angles and therefore are extremely faint and are not 
detected. In order to increase the percentage of detected 
bursts in a Monte Carlo sample of 50,000 events, I do not 
simulate bursts that lie in regions of parameter space that 


produce faint bursts; i.e. bursts with values of 9 V which 
are large compared to 6q. For each model, I plot bursts in 
the [0 o ,0 v ]-plane, showing the outline of the truncated re¬ 
gion. The models in Lamb, Donaghy, & Graziani (2005) 
correspond to removing all bursts above the line do = 9 V . 

2.3. Comparison to Burst Data 

To test the viability of various jet models, I compare 
my results against several available datasets. Consider¬ 
ing the sample of bursts with known redshifts (localized 
by BeppoSAX , HETE-2 or other detectors) there are two 
correlations of interest between source frame quantities. 
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Amati et al. ( 2002 ) reported a correlation between 
Apeak and Eiso 

A peak = , (7) 

that has been confirmed and extended to over five orders 
of magnitude in E- lso (Sakamoto et al. 2004, 2005a; Lamb 
et al. 2005). I work with the best-fit value of Ca = 
89.1 ± 8.2 keV (fixing a = 0.5 and Ea = 10 52 ergs) from 
Lamb, Donaghy, & Graziani (2005). For each realization, 
I draw a value for Ca from a lognormal distribution with 
a width of 0.13 decades. 

Recent works (Nakar & Piran 2004; Band & Preece 
2005) have claimed that large percentages of BATSE 

bursts are incompatible with the A pea k oc A iso relation 
(but see Ghirlanda et al. (2005a), Bosnjak et al. (2005) 
and Lamb et al. (2005)). Since bursts seen away from 

the Apeak oc A iso relation may be a signature of off-jet 
events, comparison with this relation is an important test 
for these models. 

Using bursts with known redshift and well-measured 
jet-break times, Ghirlanda et al. (2004) have found a 
second relation 

Apeak = Cg > ( 8 ) 

with current best-fit values of f3 = 0.69 ± 0.04, Cq = 
250 ± 100 keV and Eq = 3.8 x 10 50 ergs (Ghirlanda et 
al. 2005b). In §3.3 I describe the observational implica- 

tions of models that satisfy both the A pea k oc A is ' 0 and 
Apeak oc A^ relations. 

I define the value of A 7 measured using the method of 
Frail et al. (2001) to be A™ f = Aj so (l — cosier), where 
$br = max(0 v ,0 o )- The presence of off-jet emission (or 
any non-uniform profile) implies A™ f ^ A* rue ; in fact, 
A“ f varies with 9 V for a given jet, while A* rue is an in¬ 
trinsic property of the jet. 

The burst properties presented in the source frame 
(Ai so or A' nf against A pea k, or cumulative distributions 
of these quantities) are essentially those presented by 
Ghirlanda et al. (2004), augmented where appropriate 
with events compiled by Friedman & Bloom (2004), more 
recent fits to HETE-2 data from Sakamoto et al. (2005a), 
and a few recent Swift bursts with fits reported in GCN 
Circulars (Golenetskii et al. 2005a,b,c). 

I also consider the larger sample of HETE-2 localized 
bursts with and without known redshifts (Barraud et al. 
2003; Sakamoto et al. 2005a). This sample also shows 
a prominent hardness-intensity correlation, although it 
is broader than the source-frame correlation. This sam¬ 
ple has the advantage of having many more XRFs than 
the sample with known redshifts. Figure 3 shows both 
the observer frame and source frame datasets, with the 
relevant correlations. 

3. RESULTS 

Here I explore the relative importance of off-jet rela¬ 
tivistic kinematics for six top-hat variable opening-angle 
jet models. In what follows, bursts that are detected 
by the WXM for which 9 V < 9 0 are shown as blue dots 
(on-jet events), while detected bursts for which 9 V > 9 0 
are shown as green dots (off-jet events). Bursts in the 
simulation which are not detected are shown as red dots. 




Fig. 3.— Top panel shows the sample of bursts localized by 
the WXM on HETE-2 in the [Sg(2-400 keV),-E°^“ k ]-plane, broken 
down into XRFs, X-ray-rich GRBs and GRBs (data taken from 
Sakamoto et al. (2005a) and Barraud et al. (2003)). Bottom panel 

shows the E pea k oc Tv)'/ and E pea k oc e) relations in the source 
frame (data taken from Ghirlanda et al. (2004), Friedman & Bloom 
(2004), Sakamoto et al. (2005a) and recent Swift bursts). 


3.1. Yamazaki, I oka & Nakamura (2004a) Model 

In the first model I consider (Y04), I adopt the param¬ 
eters from Yamazaki, Ioka & Nakamura (2004a). This 
model attempts to explain classical GRBs in terms of 
the variation of jet opening-angles, while XRFs are inter¬ 
preted as classical GRBs viewed off-jet. The Y04 model 
assumes 7 = 100 and draws 9q values from a power-law 
distribution given by fo d9o oc 9q 2 dOg, defined from 0.3 
to 0.03 rad. A* rue is drawn from a narrow lognormal 
distribution centered on 1.2 x 10 51 erg. 

The upper-right panel of Figure 4 shows that the stan¬ 
dard A 7 value and a narrow range of opening angles is 
sufficient to explain the population of classical GRBs. By 

1 /2 

construction, the on-jet events follow the A pea k oc A is ' 0 
relation. Off-jet emission from similar bursts viewed at 
much larger 9 V accounts for the population of green off-jet 
points that lie above the A peak oc A is ' 0 relation. Off-jet 
events account for 34% of the total detected bursts (see 
Table 1). For different values of 9 V , these bursts gener¬ 
ally move along trajectories in the [Ai so ,A pea k]-plane that 
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Fig. 4. — Model Y04, which uses the parameters of Yamazaki, Ioka & Nakamura (2004a), where GRBs are seen on-jet and XRFs 
are explained by off-jet relativistic kinematics. Bursts detected on-jet (blue), off-jet (green) and non-detected bursts (red) are shown in 
the [0o,0v]-plane (left), [F , i so ,-F , pea k]-plane (center) and the 5s(2-400 keV)]-plane (right), with observed data from BeppoSAX and 

HETE-2 overplotted. 


1/2 

have a flatter slope than the Speak oc Ej a relation (see 
Figure 2). The observed off-jet bursts (green points) are 
consistent with only a few observed bursts, and for the 
most part, represent a population of events not seen by 
current instruments. 

In particular, the middle and right panels of Fig¬ 
ure 4 show that the HETE-2 XRFs are not easily ex¬ 
plained as classical GRBs viewed off-jet, as this model 
posits. The two XRFs with known redshifts lie along 

the .Speak oc F'iso relation, and furthermore, the larger 
sample of HETE-2 XRFs without known redshifts do not 
fall in the region of the ,Se] -plane expected for this 

model (they lie at lower, rather than higher, S°g| k values 
for a given Se). 

3.2. Top-Hat Variable Opening-Angle Jet Models 

The next two models seek to explain both GRBs and 
XRFs by a wide distribution of jet opening-angles (see 
Lamb, Donaghy, & Graziani (2005) for more details and 
discussion). These models generate XRFs that obey the 

-Speak oc E-J a relation by extending the range of possible 
jet opening solid-angles to cover five orders of magni- 
tude. Hence, XRFs that obey the S pea k oc E is ' 0 rela¬ 
tion are bursts that are seen on-jet, but have larger jet 
opening-angles. Here I add the presence of off-jet rel¬ 
ativistic kinematics to this picture. These models gen¬ 
erate a significant population of off-jet events, although 
increasing 7 reduces the fraction of off-jet bursts in the 
observed sample. 

I draw Ho = 27t(1 — cosf?o) values from a power-law 
distribution given by fo dtto oc H(( 2 dtt 0 , defined from 27t 
to 27 t x 10 -5 sr. S* rue is drawn from a narrow lognormal 
distribution centered on 1.2 x 10 49 erg. The lower cen¬ 
tral point for the S 4rue distribution is a requirement for 
including in a unified model those events with measured 
E iso values that are smaller than the usual standard en¬ 
ergy of ~ 10 51 ergs. I consider 7 = 100 (THVOAl, 
Figure 5) and 7 = 300 (THVOA2, Figure 6 ). 

As can be seen from the upper-left panels of Figures 5 
and 6 , the relative importance of off-jet events increases 
for models with a population of very small opening- 


angles. This is mainly due to the fact that for con¬ 
stant S 7 , narrower jets will have larger Sj so values. Such 
bursts will be brighter and therefore detectable at larger 
6 V . Narrow jets are also more likely to be viewed off-jet 
than wider jets. Finally, smaller values of 9q give rise 
to trajectories in the [Si so ,S pea k]-plane that differ more 

conspicuously from the Speak oc EJ 0 relation. 

The large population of very small jet opening angles 
in these models also implies a much larger true burst 
rate than in model Y04. As Table 1 shows, for these 
two models the WXM will detect 1 out of 2570 and 6470 
bursts respectively, as compared to 1 out of 150 bursts 
for the Y04 model . 2 

For larger values of 7 , the emission curves in Figure 1 
drop off faster away from the edge of the jet. Comparing 
the upper-right panels of Figures 5 and 6 illustrates the 
fact that larger values of 7 reduce the percentage of off- 
jet events observed by the WXM, and consequently the 

1 /2 

percentage of bursts seen away from the S pea k oc E is ' 0 
relation. This population of events is fairly conspicuous 
in the 7 = 100 case (83% of detected bursts are off-jet 
events) and less so for 7 = 300 (56%). 

Furthermore, this figure highlights the fact that ob¬ 
served XRFs are more easily explained by wide opening- 
angle jets than by off-jet emission. No matter the value 
of 7 , observed XRFs inhabit different regions of the data 
planes than do off-jet events. 

3.3. Models Matching the S pea k oc E@ Relation 

Reproducing the Speak oc E@ relation (Ghirlanda et al. 
2004) is an important test of any jet model. Discovery 
of this relation post-dated Lamb, Donaghy, & Graziani 
(2005), and so was not addressed in that paper. Here I 
construct models that satisfy this relation, in addition to 
the constraints of the earlier models (for more details see 
Donaghy, Lamb, & Graziani (2005c)). 

The Speak oc sV 2 and S pea k oc E@ relations can be 
mutually satisfied in the top-hat variable opening-angle 

2 In calculating the detected fractions, I correct for the trunca¬ 
tion in the simulation parameter space described in §2. See §4 for 
more details. 
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Fig. 5.— THVOA1 model, which explains both GRBs and XRFs by a wide distribution of jet opening-angles, using 7 = 100. Bursts 
detected on-jet (blue), off-jet (green) and non-detected bursts (red) are shown in the [ 0 o, 0 v]-plane (left), [/fi so ,/f pea k|-plane (center) and 
the [F'peak’'^’B(2-400 keV)]-plane (right), with observed data from BeppoSAX and HETE-2 overplotted. 
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Fig. 6 .— THVOA2 model, which explains both GRBs and XRFs by a wide distribution of jet opening-angles, using 7 = 300. Bursts 
detected on-jet (blue), off-jet (green) and non-detected bursts (red) are shown in the [ 0 o, 0 v ]-plane (left), [Fj so ,/f pea k]-plane (center) and 
the [F^s k ,5 £ ( 2 -400 keV)]-plane (right), with observed data from BeppoSAX and HETE-2 overplotted. 


jet model by imposing a relation between flo and E 7 . 
The following expressions hold only for the on-jet events; 
the off-jet events behave differently. Combining Equa¬ 
tions 7 and 8 with the definition 


EL 


E. 


7 


I find 


E pea k — CA 


Slo/27 


E-y ■ 2n 
Ea • flo 


1 — cos t 




(9) 


( 10 ) 


and solving for flo gives the following expression 


^0 = 2 tt (1 - cos 6q) = 2?r f 


1/a 


e(° 

7 


-0)/a 


( 11 ) 

I consider two different models in which flo and E 7 rue 
are specified by imposing this relationship. The natural 
minimum value for the distribution of E 7 rue is the point 


E* 


at which E 1 = E-, 


lso , which is found to be 


E* = 


( C a Eq\ 

\Cc,E<a) 


( 12 ) 


Current parameters for the correlations give a value of 
E* = 3 x 10 44 ergs, which is well below current observa¬ 
tional thresholds. Values of E* Tue are generated by draw¬ 
ing from a power-law distribution that gives equal num¬ 
bers of bursts per decade, and is defined from E* through 
a maximum value (E 7 rue = 3.16 x 10 51 ergs) that encom¬ 
passes the largest observed E 7 value (EJ™ f = 5.75 x 10 51 
ergs for GRB 990123). To avoid a sharp cutoff at the 
minimum and maximum values of 2 £ 7 , I include an addi¬ 
tional smearing function in the simulated value of E* lue 
with a lognormal width of 0.3 decades, flo (and hence 
do) are found via Equation 11, and the rest of the sim¬ 
ulations proceed as above. Results for 7 = 100 (GGL1) 
and 7 = 300 (GGL2) are shown in Figures 7 and 8 , re¬ 
spectively. 

Current data seems to indicate that the E P e a k oc E@ 


relation has a narrower distribution about the best-fit 
line than does the E pea k oc E is ' 0 relation (Ghirlanda et 
al. 2004). Adding a Gaussian smearing function to Ca, 
as was done in the above models, produces equal widths 

for both distributions. To broaden the E pea k oc E-J 0 re¬ 
lation, I introduce an additional Gaussian smearing func¬ 
tion into the relation between E 7 and fig- Equation 11 
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then becomes 


n 


2n 

= C^ 



El) 


1/a 

Trl a ~P)/ a 


(13) 


where Cq is centered on 1.0 and has a lognormal width 
of 0.3. Combined with the above value of the smearing in 
C A , this approach and value of Cq gives good agreement 
with both observed distribution widths. 

An immediate consequence of the flo -E 1 correlation is 
a lack of very small jet opening-angles and a narrower 
range of 9 o (compare the upper-right panels of Figures 
6 and 8 ). The bursts observed off-jet lie closer to the 

-Speak oc A is ' 0 relation than in models without the cor¬ 
relation (compare the upper-left panels of Figures 6 and 
8 ). This is due to the narrower range of 9q values men¬ 
tioned above; larger 9q values produce trajectories in the 
[i^iso^peak]-plane that closer approximate a 0.5 power- 
law. The off-jet bursts also deviate further from the 
Sp ea k oc E@ relation than from the S pea k oc SV 2 rela¬ 
tion. This is due to that fact that, given an on-jet and an 

1/2 

off-jet event that lie near each other on the S pea k oc S is ' 0 
relation, the off-jet event is more likely to have a smaller 
$br than the on-jet event, thereby giving a smaller E™ { 
value. 

In comparison with the two THVOA models, these 
models exhibit much smaller true burst rates (1 detection 
for every 165 and 197 bursts, respectively) and smaller 
fractions of off-jet events (26% and 10 % of detected 
bursts, respectively). The additional Gaussian smearing 
function, Cq, has the effect of blurring the clear sepa- 
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ration of blue and green points along the A pea k oc A iso 
relation (compare THVOAl and GGL1), but not along 
the -Speak oc E@ relation in which Cq plays no role. 


3.4. Models with flo-7 Correlations 

Finally, I investigate the effect of an additional corre¬ 
lation between fio and 7 . If narrower jets have larger 
bulk 7 values, this could reduce the importance of off-jet 
emission even further. The exact relationship between 
the bulk 7 of the material and the opening-angle of the 
jet is unknown. Extending the GGL models in the previ¬ 
ous section, I consider a simple model (GCOR) in which 
7 is given by 7 oc . I fix the normalization by setting 
7 = 1000 for the bursts with the narrowest jets. 

Figure 9 shows that imposing this correlation greatly 
reduces the percentage of bursts seen off-jet (< 7% of 
detected bursts are seen off-jet). Physically this is due 
to a combination of several effects. It is more probable to 
observe narrow jets at viewing angles slightly outside the 
jet than inside the jet, yet for such narrow jets a large 
value of 7 ensures that the detectability of such slightly- 
off-jet bursts drops off very quickly with viewing angle. 
Broader jets may have smaller values of 7 , but they are 
less likely to be observed off-jet. 

4. DISCUSSION 

Table 1 summarizes the detected fractions and off-jet 
fractions for the six models. The detected fraction has 
direct implications for the total rate of GRBs. The ratio 
of the observed rate of Type Ic supernovae to the ob¬ 
served rate of GRBs is roughly 10 5 (Lamb 1999, 2000), 


and therefore the ratio of the observed rate of Type Ic 
supernovae to the true rate of GRBs is that value times 
the detected fraction for the model. Due to their very 
narrow jets, the THVOA models have the smallest de¬ 
tected fractions of the six models presented here. For 
the model THVOA2, GRBs may comprise an apprecia¬ 
ble fraction of all observed Type Ic supernovae, but for 
all other models, the true rate of GRBs is much smaller 
than the observed supernova rate. 


Model 

Detected Fraction 

Off-Jet Fraction 

Y04 

0.00659 

= 1/151.7 

0.339 

THVOAl 

0.00039 

= 1/2571.1 

0.829 

THVOA2 

0.00015 

= 1/6467.3 

0.561 

GGL1 

0.00606 

= 1/165.0 

0.262 

GGL2 

0.00507 

= 1/197.3 

0.103 

GCOR 

0.00489 

= 1/204.4 

0.066 


TABLE 1 

Detection statistics for the six models presented in this 
paper. The second column shows the true rate of 
detection for all bursts in the model. Note that this is 

NOT THE SAME AS THE PERCENTAGE OF BURSTS DETECTED IN THE 
SAMPLES OF 50, 000 DUE TO THE TRUNCATION DESCRIBED ABOVE. I 
CORRECT FOR THAT TRUNCATION IN CALCULATING THE DETECTION 
FRACTION. THE THIRD COLUMN SHOWS THE FRACTION OF ALL 
DETECTED BURSTS THAT ARE SEEN OFF-JET. 


In the figures above, I compare theoretical models em¬ 
ploying the WXM detector threshold with data compiled 
from many instruments with varying detector sensitivi¬ 
ties. To assess the effect of the detector thresholds in the 
simulations, I compare the predicted distributions of ob¬ 
served bursts for the two most successful models (GGL2 
and GCOR) using the WXM instrument and using the 
GRBM instrument on BeppoSAX. For these two models, 
I compare the predicted distribution of observed bursts 
employing a given instrument threshold only with the 
data from that instrument. 

Figure 10 shows that the higher triggering threshold 
for the GRBM instrument on BeppoSAX prevented that 
mission from promptly localizing the fainter, low-A pea k 
XRFs. In contrast, rapid HETE-2 localizations of XRFs 
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have provided evidence that the Apeak oc A is ' 0 relation 
extends to lower A pea k values, but HETE-2 has detected 
fewer high-A pea k bursts than BeppoSAX. Therefore, the 
models presented in this paper employing the WXM 
threshold are a good match for the full-range of observed 
GRB characteristics. 

Figures 11, 12 and 13 summarize the six models by 
comparing the observed data sets against the model cu¬ 
mulative distributions of E iso , A peak , A" lf , 9 hr , A£^ k , 
and Se (2-400 keV). Although a great deal of information 
is lost in projecting the 2 -dimensional distributions from 
the above figures onto each axis separately, some useful 
information can still be obtained from these curves. 

Figure 11 shows the cumulative distributions for Ai so 
and Ap ea k, again separately comparing models using ei¬ 
ther the WXM or GRBM instrumental threshold with 
data obtained from that instrument. Due to the small 
size of the current datasets, the THVOA, GGL and 
GCOR models would not be judged as inconsistent with 
the data by the KS test. Model Y04 is here seen to over¬ 
produce brighter, high-A pea k events, and under-produce 
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Fig. 7.— GGL1 model, which picks parameters to match both the i£ pea k oc and i£ pea k oc E^ relations, using 7 = 100. Bursts 

detected on-jet (blue), off-jet (green) and non-detected bursts (red) are shown in the [#O) 0 v]-plane (left), [i£,i£ pea ]J-plane (center) and the 
[-^peak’^-E(2-400 keV)]-plane (right), with observed data from BeppoSAX and HETE-2 overplotted. 
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Fig. 8 .— GGL2 model, which picks parameters to match both the Speak 77 and Speak oc S^ relations, using 7 = 300. Bursts 
detected on-jet (blue), off-jet (green) and non-detected bursts (red) are shown in the [Vj/i v ]-pla.ne (left), [S,S pea k]-plane (center) and the 
[S°k°k,‘S , £(2-400 keV)]-plane (right), with observed data from BeppoSAX and HETE-2 overplotted. 



low-Sp eak XRFs. 

Figure 12 shows the cumulative distributions for 
and (?br 1 again separately comparing models using either 
the WXM or GRBM detector threshold with data ob¬ 
tained from that detector. This dataset is even sparser 
than that of Figure 11, but the figure does highlight 
the effect of rescaling the standard energy downward to 
10 49 ergs for the two THVOA models. The GGL and 
GCOR models avoid this problem by incorporating the 
-Epeak oc E@ relation, and thereby avoid the large dis¬ 
parity with the and 6*br distributions seen for the 
THVOA models. 

Figure 13 shows the cumulative distributions for E° 7 k 
and Se (2-400 keV), comparing models using the WXM 
detector threshold with the larger dataset of all bursts 
localized by HETE-2. Again, the THVOA, GGL and 
GCOR models are all reasonably consistent with the 
data. There is some hint that the GGL and GCOR mod¬ 
els produce more bright, high-Speak bursts than were 
seen by HETE-2, but this may be the result of trying 
to match both the BeppoSAX and HETE-2 populations 
with one model. 


To summarize, models GGL2 and GCOR are the most 
successful at matching the 2-dimensional distributions 
discussed above as well as the individual cumulative dis¬ 
tributions. Model Y04, which seeks to explain XRFs as 
off-jet GRBs, is unable to match the distributions of ob¬ 
served XRFs. The THVOA models produce large num- 
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bers of detectable bursts away from the A peak oc EJ 0 
relation which are not seen in current datasets, and re¬ 
quire a re-scaled standard energy of 10 49 ergs, which is 
inconsistent with current afterglow theories (Panaitescu 
& Kumar 2001). 

5. CONCLUSIONS 

Bursts with known redshifts have been found to obey 

the K P e a k oc S iso relation, and a large population of 
bursts away from this relation is not readily apparent 
in current datasets. In particular, the limited sample of 
XRFs with known redshift information is consistent with 
an extension of the S peak oc E-J 0 relation to over 5 orders 
of magnitude in Sj so . Liang, Dai & Wu (2004) found that 

the l? pea k oc EJ o relation holds internally within a large 
sample of bright BATSE bursts without known redshift, 
perhaps indicating that the relation is a signature of the 
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Fig. 9.— Model GCOR where 7 oc and 7 = 1000 at the minimum value of Bursts detected on-jet (blue), off-jet (green) and 
non-detected bursts (red) are shown in the [0o 5 0v]-plane (left), [i?,.Ep ea k]-plane (center) and the [F , p^ k ,S , £:(2-400 keV)]-plane (right), with 
observed data from BeppoSAX and HETE-2 overplotted. 




Fig. 10. — The effect of different detector thresholds for two models, GGL2 (upper row) and GCOR (lower row). Bursts are shown in 
the [F 1 ,Speak]-plane, comparing models employing the WXM threshold against WXM data points (left column) and models employing the 
BeppoSAX threshold against BeppoSAX data points (right column). 


physics of the emission mechanism. However, recently 
some authors have argued that the Speak oc S iso rela¬ 
tion may be an artifact of some unknown selection effect 
arising from the process of determining the burst red- 
shift, and that 25% (Nakar & Piran 2004) to 88 % (Band 
& Preece 2005) of the BATSE bursts may be inconsis- 

1 /2 

tent with the S pea k oc S iso relation. However, these 
results are controversial (Ghirlanda et al. 2005a; Bosn- 


jak et al. 2005; Lamb et al. 2005) and depend sensitively 
on the quality of the spectral fit that generates the S°g® k 
parameter. I therefore regard the question of the per¬ 
centage of BATSE bursts that are inconsistent with the 
Speak oc S iso relation to be an open one. 

For low values of 7 , top-hat variable opening-angle jet 

models predict a sizable population of detectable, off-jet 

1 /2 

bursts that lie away from the S pea k oc S iso relation and 
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Fig. 11.— The cumulative distribution in Fj so (left column) and Speak (right column) for detected bursts in the six models explored 
in this paper. The upper row compares models employing the GRBM detector with data points observed by BeppoSAX.. The lower row 
compares models employing the WXM detector with data points observed by HETE-2. The data used for the S pea k and E lso distributions 
are the same as in Figure 3b, broken down by detector. 


that are not seen in current data sets. It may be that 
such bursts are in fact present in the BATSE catalog 
and form a population of bursts that do not obey the 
-Speak oc E.J o relation. On the other hand, if such a 
population is found not to exist, it implies that the bulk 7 
of the jet is large (~ 300 — 1000). For models that include 
the Speak oc E@ relation, the off-jet burst population lies 

closer to the S pea k oc S is ' 0 relation than for other models, 
but a similarly discrepant population of off-jet bursts is 
found to lie above the S pea k oc E@ relation, leading to 
similar conclusions. 

Regardless of the size of the population of off-jet bursts, 
it seems unlikely that such a population can make up the 
bulk of the XRFs observed by HETE-2. The larger sam¬ 
ple of HETE-2 bursts without redshift information con¬ 
tains XRFs which lie toward smaller E°^ k values than is 
predicted by the off-jet emission model and hence are not 
easily explained as classical GRBs viewed off-jet. Mod¬ 
els in which XRFs are the product of larger jet opening- 
angles seem to better match the observed distributions 
of XRF properties. This seems to match the evidence 
arising from X-ray afterglows of XRFs. Granot et al. 
(2005) calculated the afterglow light curves predicted by 
various models of burst emission seen off-jet. They find 
that a general feature of off-jet afterglows is an initial 
rising light-curve that peaks at about the jet break time 
and then declines rapidly, similar to an on-jet event. Af¬ 


terglows with initially rising components have not been 
observed. In particular, XRFs with well-observed X-ray 
afterglows, for example XRF 020427 (Amati et al. 2004) 
and XRF 050215b (Sakamoto et al. 2005b), have after¬ 
glow light curves that join smoothly onto the end of the 
prompt emission and that show no evidence of a jet break 
for many days after the burst, implying large jet opening- 
angles. 

It is straightforward to arrange for top-hat variable 
opening-angle jet models to match the empirical .Speak oc 
EP relation, and such models also provide a natural ex¬ 
planation for XRFs. Figures 7 and 8 illustrate the conse¬ 
quences of adopting the correlation between and E* lue 
that ensures that on-jet events obey the S pea k oc E@ re¬ 
lation. Most importantly, incorporating the E pea ^ oc 
relation in the THVOA model removes two of the main 
drawbacks of the THVOA model presented in Lamb, 
Donaghy, & Graziani (2005). The requirement of very 
small (~ 2 °) jet opening angles to explain the largest 
Ei so values was criticized (Stern 2003) as being difficult 
to achieve in a hydrodynamic jet. The high end of the 
Sj so distribution is here explained by jets with moderate 
opening angles but larger E-, values. The need to re-scale 
the central value of E* lue downward to ~ 10 49 ergs to in¬ 
corporate the XRFs in a unified model was criticized as 
being difficult to reconcile with afterglow models. The 
-Epeak oc E@ relation naturally produces a range of E 1 
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Fig. 12.— The cumulative distribution in E (left column) and #br (right column) for detected bursts in the six models explored in this 
paper. The upper row compares models employing the GRBM detector with data points observed by BeppoSAX. The lower row compares 
models employing the WXM detector with data points observed by HETE-2. The data used for the E ^ nf distribution are the same as in 
Figure 3b, broken down by detector. The data used in the $br distribution is the dataset from Ghirlanda et al. (2004). 




Fig. 13.— The cumulative distribution in (left) and S#(2-400 keV) (right) for detected bursts in the six models explored in this 

paper. These panels compare models employing the WXM detector with the larger dataset of all bursts localized by HETE-2. The data 
used for these distributions are taken from Sakamoto et al. (2005a) and Barraud et al. (2003). 


values that extends down into the XRF regime. 

Matching the U pea k oc E'?, relation also mitigates the 
problem of a large population of bursts seen away from 

the Speak oc EJ 0 relation in the low 7 case. In these 
models the off-jet events hew more closely to the S pea k oc 
1 /2 

E-J Q ~ relation, but a similar problem arises in that these 
off-jet events are seen away from the S pea k oc E@ relation 


instead. A possible way out of this dilemma might be to 
impose a relationship whereby narrower jets have larger 
bulk 7 values and broader jets have smaller 7 values. 
Using a simple model where 7 oc 1/f2o results in a sub¬ 
stantial reduction in the percentage of detected bursts 
seen off-jet. 

Yamazaki, Ioka & Nakamura (2004b) have proposed a 
multiple subjet model for unifying short and long GRBs, 
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X-ray-rich GRBs and XRFs. The model employs emis¬ 
sion from multiple subjets (seen off-subjet) to explain 
X-ray-rich GRBs and XRFs. The authors performed 
Monte Carlo simulations for a universal multiple subjet 
model and find that the results are consistent with the 
-Epeak oc E-J a relation, albeit with considerable scatter 
(see Figure 4 in Toma et al. (2005)). 

There are two reasons why the multiple subjet model 
better satisfies the .E pea k oc E-J 0 relation. First, the 
authors choose 7 = 300, which satisfies the constraint 
on 7 that I find above. The behavior of bursts viewed 
outside the envelope of subjets should approximate the 
top-hat models considered in this work. If the authors 
had adopted a lower value of 7 , the model would have 
produced a large number of bursts that lie away from the 
E pe ak oc EE q relation; i.e., they would have encountered 
the same problems as those of model THVOA1. Second, 
since 7 = 300, the spectrum for each line of sight is 
dominated by that of the closest subjet, and since there 
are many subjets, each line of sight lies very close to 
at least one subjet, mitigating the effects of relativistic 
kinematics produced by viewing a subjet well off the jet. 

Finally, Toma et al. (2005) find that for values of 
7 lower than 300, the ratio between GRBs, X-ray-rich 
GRBs and XRFs becomes highly skewed toward hard 
GRBs, in contradiction with the HETE-2 results. Thus, 
all of the results in Toma et al. (2005) support the re¬ 
quirement I find in this paper that large gamma values 
are needed in order to match the observed data for XRFs, 
X-ray-rich GRBs, and GRBs. 

Off-jet relativistic kinematics will be important in non- 
uniform jets as well as in top-hat jets. Models employing 
Gaussian (Zhang et al. 2004) or Fisher-shaped (Donaghy, 
Lamb, & Graziani 2005a) jets rely on the exponential fall 
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off of the emissivity with viewing angle to match the wide 
spread of observed burst quantities. By including off-jet 
emission in these models, the exponential fall off will be 
dominated at some angle by the power-law fall off due 
to relativistic kinematics, thereby broadening the emis¬ 
sivity distribution and reducing the range of generated 
Ei so values. Graziani et al. (2005) showed that different 
underlying burst profiles may have radically different ob¬ 
servational distributions. We hope to use population syn¬ 
thesis Monte Carlo simulations to further explore these 
models in future work. 

In conclusion, off-jet emission from collimated GRB 
outflows should exist simply as a consequence of rela¬ 
tivistic kinematics. Monte Carlo population synthesis 
simulations of top-hat shaped variable opening-angle jet 
models predict a large population of off-jet bursts that 
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are observable and that lie away from the E pea k oc EE 0 
and Apeak oc E@ relations. Such off-jet events are not ap¬ 
parent in current datasets. These discrepancies can be 
removed if 7 > 300 for all bursts or if there is a strong 
inverse correlation between 7 and Do- The simulations 
show that XRFs seen by HETE-2 and BeppoSAX cannot 
be easily explained as classical GRBs viewed off-jet, and 
are more naturally explained as jets with large opening- 
angles. 
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